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We formulate non-relativistic classical and quantum mechanics in the non-commutative two dimen¬ 
sional plane. The approach we use is based on the Galilei group, where the non-commutativity is 
seen as a central extension upon identification of the boost generators with the position operator. 
We perform a systematic study of the free particle, defined by the symmetries of the space time, 
which include the no-commutativity. The symmetries at the classical level are analyzed in terms of 
Noether’s theorem. Canonical quantization is presented and the representation of the corresponding 
Heisenberg algebra is obtained. The path integral representation and Wigner distribution function 
in phase space are also discussed. We work out, both at the classical and at the quantum level, the 
harmonic oscillator, avoiding the use of the conventional non-canonical transformation that leads to 
a momentum dependent potential. We use Einsteins’ model for a solid to corroborate that, accord¬ 
ing with intuition, as a consequence of the space fuzziness the entropy is a growing function of 0 in 
the low temperature regime. 
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I. INTRODUCTION 

Non-commutativity (NC)^ can be interpreted either, as 
an intrinsic property of space time Q , or as a low-energy 
consequence of the fundamental string theory Q . In both 
cases. Non commutative quantum mechanics (NCQM) is 
seen as a laboratory where the properties of NC field 
theories can be studied. This has however lead to am¬ 
biguities of what is meant by NC mechanics, in fact dif¬ 
ferent points of view exist in the literature which seem 
to be equally relevant. Among the more frequently used 
are the one based on the Moyal product, i.e. assume 
that the dynamics is described by the conventional equa¬ 
tions but using everywhere Moyal products, and the sec¬ 
ond approach where one suppose that the only way NC 
enters is through the Poisson Brackets (commutators in 
QM) of the position operators. How is the NC related to 
the space-time, i.e. can we identify the non-commuting 
operators with the position? and if so, is the resulting 
theory consistent, do we have precise rules to describe 
the physics in different reference frames?. From our view 
point one should clearly state the assumptions made to 
define the NC mechanics we use, in particular, one should 
know what is the effect of NC on a free particle. 

In this paper we adopt the first view, namely we assume 
NC is intrinsic property of space time and require that 
NC be incorporated in a consistent way. Indeed, an im¬ 
portant property of any formulation of the mechanics is 
the symmetry of the underlying space-time, i.e. the rel¬ 
ativity group on which it is based. In fact, we will fol¬ 
low the procedure used to formulate field theory where 


the free particle is treated using the symmetries of the 
space time and then, once a consistent framework is avail¬ 
able, interactions are introduced using further arguments 
(renormalizability, gauge invariance, SUSY). Thus, we 
first consider in detail the free particle taking as basis a 
Hamiltonian description based on a symplectic structure 
- determined by the symmetries of the NC space - which 
we assume is the Galilei group. It is important to remark 
that already at this level we have some constraints since 
consistency of NCM with Galilei group requires: 

• to work in the two dimensions, 

• although {xi,Xj} ^ 0 necessarily {pi,pj} = 0. 

An important point to remark is that in our approach, 
the Poisson Brackets (PB) (or commutators in the quan¬ 
tum theory) of the coordinates is always non zero, we 
avoid any non-canonical transformation leading to van¬ 
ishing PB, otherwise it would be inconsistent with the 
assumed Galilei algebra. Thus, for us, the interaction 
of a charged particle in two dimensions with a perpen¬ 
dicular constant magnetic field is not a typical example 
of NC since in such a case the free particle reduces to 
the standard commutative problem which can be consis¬ 
tently quantized and, after that, the interaction can be 
introduced through minimal coupling. We will comment 
on these points below, but we first refer to some of the 
existing work in the literature. Since there are so many 
publications in this area, we will restraint to those of 
direct interest to our approach, and still, apologize for 
undeliberate omissions. 


‘Electronic address: vaquera@fisica.ugto.mx, lucio@fisica.ugto.mx 
^ In the following we will use NC to stand for non-commutative and 
M, C and Q for mechanics classical and quantum respectively. 


From our point of view the formulation of NCM that 
not only is consistent with, but actually is based on, the 
Galilei algebra is of particular interest. In such an ap¬ 
proach - which is valid only in two dimensions |3) 0 “ 
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the NC appears, together with the mass, as two cen¬ 
tral extensions of the algebra. This was recognized by 
authors in Q who build a model that provides a realiza¬ 
tion of the symmetry. They also discuss the quantization 
of the model using the formalism of constraints and a 
lagrangian including higher derivatives. In spite of the 
completeness of the work in Q, there are several top¬ 
ics that deserve further discussion. Thus, in the present 
paper we make a summary of non-commutative classical 
mechanics NCCM in two dimensions. Besides the Hamil¬ 
tonian formalism, the first order Lagrangian formulation 
of NCCM and the relation between the Hamiltonian and 
Lagrangian are also discussed following the approach in 
0 . The equivalence between those formalism natu¬ 
rally leads to the conclusion that two dimensional NCCM 
can always be treated as a system with second class con¬ 
straints, at least if interaction with gauge fields are not 
introduced 0- We discuss Noether’s theorem and start¬ 
ing from the symmetries we derive the Hamiltonian for 
the free particle on the NC plane. Although this appears 
to be a trivial and long exercise, in fact it shed some light 
on some of the questions previously formulated. Indeed, 
besides helping in the identification of the Hamiltonian 
(there is not consensus on this point, for example accord¬ 
ing to the approach in 0 the dynamics of a free particle 
on a NC space is equivalent to the study of the dynamics 
of this particle with charge q on the usual commutative 
space in the presence of a magnetic field), the symme¬ 
tries determines the Hamiltonian as the time evolution 
operator, providing thus some support to the whole ap¬ 
proach. All along the the paper we emphasize the role 
played by symmetries, this is relevant since the proper 
formulation and solution of some of the problems in QM 
can be traced to the appropriated understanding of the 
classical symmetries. 

With the formulation of the mechanics previously dis¬ 
cussed the passage from classical to quantum mechanics 
is straightforward. At the quantum level, using canon¬ 
ical quantization, we analyze the representation of the 
Heisenberg algebra as a way to implement the canonical 
quantization. As explained below, the idea is to imple¬ 
ment representations for the operators which are consis¬ 
tent with the commutation relations, but also determine 
what is the relevance of the commutation relation for the 
states. The advantage of this approach is that, based on 
the commutation relations, one can implement and com¬ 
pare different representations. An example where a judi¬ 
cious choice of the base is exploited is seen in 0 where 
the periodicity of the oscillator in the non-commutative 
plane is analyzed, although neither the spectrum nor the 
eigenfunctions of the complete set of observables are de¬ 
rived. Previous work include also an analysis of possible 
realizations of the operators, which is a related but not 
equivalent problem |in|. . We derive expression for 

non-equivalent representations of the Heisenberg algebra 
extending previous analysis to NCQM ^3- The repre¬ 
sentations include differential representations of the op¬ 


erators (position and momenta) in different basis as well 
as gauge fields that follow from the structure of the alge¬ 
bra and that are relevant in the description of non trivial 
manifolds. As far as we know, ours is the first time a 
detailed analysis of the representations of the Heisen berg 
algebra is presented in the two dimensional NC case 

An alternative procedure to study the quantum proper¬ 
ties of a system is through the phase space Wigner dis¬ 
tribution function. Besides the intrinsic interest on this 
formulation, the use of non-commutativity in the field of 
quantum optics 0 , where Wigner distribution function 
is a common tool, is a further motivation for its general¬ 
ization to the non-commutative case. We are not aware of 
exhaustive work along this line, the difference with exist¬ 
ing literature being again the approach 0 We provide 
a compact exact expression for the Wigner distribution 
function, again in this case we do not need to consider an 
expansion in 0, the NC parameter. For completeness we 
include a brief discussion on the path integral formula¬ 
tion, where the results of the canonical quantization are 
used, in particular the wave functions that permit to con¬ 
nect different basis. Our results are analogous to those 
obtained in 0 

As examples of the application of the formalism, we work 
out in detail two problems, both at the classical and at 
the quantum level: the free particle and the harmonic 
oscillator. The free particle can be considered the phys¬ 
ical system where the symmetries of the space time are 
realized, for that reason we characterize the free particle 
with the Hamiltonian that is compatible with the corre¬ 
sponding Galilei group. On the other hand, several pa¬ 
pers [hi HTl have dealt with the harmonic oscillator 
at the quantum level. The case of the two dimensional 
harmonic oscillator in the presence of a constant mag¬ 
netic field has been considered 00 and it has been 
remarked that two phases exist, in one of these (B ^ 0) 
the symmetry group is SU{2) while in the other (B = 0) 
the symmetry group is 5*17 (1,1). In our approach we 
show in CM, that although SU(2) is a symmetry of the 
commutative harmonic oscillator, this is not true in the 
NC case. We analyze the surviving symmetry and, at the 
quantum level, determine the eigenfunctions common to 
the Hamiltonian and angular momentum. 

Before concluding this introduction in this paragraph we 
make a summary of our approach. We consider the me¬ 
chanics, both classical and quantum, for a free particle 
imposing consistency with the symmetry of the space 
time, which is assumed to be Galilei group including two 
central extensions (the mass and the parameter of NC) 

. We completely avoid the use of non canonical trans¬ 
formations and assume that the only way NC enters is 
through the commutation relations of the position oper¬ 
ators plus those that the requirement of consistency de¬ 
mands. A potential V{x)^ independent of the NC param¬ 
eter 0, is added to the Hamiltonian of the free particle. 
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which has been previously determined through the sym¬ 
metry. In QM we assume the validity of the Schrodinger 
equation and use the time translation generator, i.t. the 
Hamiltonian. The differential operator associated to the 
Hamiltonian is obtained from the classical one plus any 
of the four representations of the position and momenta 
consistent with the Heisenberg algebra. We have tried 
to make the paper as self contained as possible, avoid¬ 
ing however unnecessary details. We organized the main 
body of the manuscript in six parts. Section two is de¬ 
voted to classical mechanics, the third to the canonical 
quantization. Sections four and five deal with alternative 
quantization procedures. Examples in classical mechan¬ 
ics are discussed in section two, in QM in section six and 
we end with a summary of the contributions of this work. 


II. CLASSICAL MECHANICS 

We start with a brief summary of classical mechanics. 
This will allow us to introduce the notation and also to 
clarify the role played by symmetries in the formulation 
of NC mechanics. In the Hamiltonian formalism the de¬ 
scription of a system with n degrees of freedom is deter¬ 
mined by 2n first order differential equations involving 
2n independent variables that, here and thereof, we de¬ 
note by Za (a = 1,2, ...2n). The symplectic structure J 
associated to this formalism defines the Poisson brackets 




dz>^ ’ 


( 1 ) 


where summation over repeated indices is understood. 
The parenthesis are real, antisymmetric and linear. Fur¬ 
thermore they satisfy the Leibiniz rule and the Jacobi 
identity: 

• {f,9h} = {f,9}h + 9{f.h} 


• {f,{9,h}} = {{f,g},h} + {g,{f,h}}. 

According to o the entries of J are: 

= = ( 2 ) 

For later reference we introduce the inverse of the sym¬ 
plectic structure J by means of the relation: 

= (3) 


the conjugated momenta with the last n {pi = Zi,i = 
n-|-1, n + 2, ...2n), then the information on the symplectic 
structure J is contained in the Poisson brackets: 


{9a^ 9 / 3 } — ^a/3j {QouP/s} — ^a)3j {PonPf^} — ©ct/Jj (5) 


So far = Oapig^p) is an antisymmetric field = 
—Ofja) which, in general, may depend upon the position 
and the momenta. However, if we set 0a/3 = 0, the Jacobi 
identity (applied to a momentum pi and two coordinates 
qj,qk), implies the relation: 


9% 

dqk 


= 0 , 


( 6 ) 


therefore Oap = 0c</3(p)- At this point it should be clear 
that, by allowing a general enough symplectic structure, 
NCM can be described within this formalism. 


On the other hand, the lagrangian formulation is rele¬ 
vant in discussing the symmetries. For that reason we 
introduce the first order Lagrangian Q: 

L{z,z,t) = Ka{z)z°‘- H{z,t), (a=l,---,2n) (7) 


where Ka is a vector potential in phase space for the 
inverse of the symplectic structure, i.e.: 


■^CX.13 


dKp dK^ 
dz°‘ ~ dz0 


( 8 ) 


in order to show the equivalence of this Lagrangian to the 
Hamiltonian formalism previously introduced, we apply 
the variational principle to the action associated to this 
Lagrangian: 


=sj^' {K^{z)z^ - H{z, t)) 


dt 


( 9 ) 




rtf r 


iU/SaZ 


dH 


Sz^dt + {Kc,Sz‘^)\lt. 


Therefore, if the variation at the end points are such that: 


KJz‘^\t^=K^Sz%^, 


the variation of the action Eq. 0 implies: 


W/3ai“ 


dH 
dz!^ ’ 


( 10 ) 


( 11 ) 


The equations of motion for a system described by the 
Hamiltonian H{z), in a phase space with symplectic 
structure J, are given by: 

i“ = K,H(0)} = a,/I = l,...,2n. (4) 

If one identifies the generalized coordinates with the first 
n variables {qi = Zi,i = l,2,...n) of phase space and 


which are nothing but the Hamilton equations of motion 

Eg.CT: 

i“ = {z“,H(z)} = (12) 

If uj(z) has constant entries, Ka can be written as: 

Ka = ^Z^UJfSa- (13) 
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and, under these conditions, the Lagrangian reduces to: 

(14) 


L{z,z,t) = ^Z^UJ/SaZ^^ - H{z,t), 


With the lagrangian formulation at hand, Noether’s the¬ 
orem is formulated in the conventional way, for complete¬ 
ness we quote the result 

If under a group of order r, of continuous transformations 
t ^ t' = > z'“(z, t), the action is invariant up 

to surface terms i.e. 


S[z’] = S[z] + 


dA{z) 

dt 


dt. 


(15) 


then, for every classical solution to the equations of mo¬ 
tion, r functions 7 a (A = 1, 2, ...r) of the dynamical vari¬ 
ables are conserved: 


7a = -Xx{H- z^^K^) - Aa. 


(16) 


the quantities (j), x and A are given by the infinitesimal 
transformations which are a symmetry of the action: 


St = e^x\{t), 

5A = e^Ax, 


(17) 


where e"^ (A = 1, • • • , r) are constant, infinitesimal pa¬ 
rameters each of which is associated to a given transfor¬ 
mation. It proofs convenient to introduce the notation 
Q = e^7A (no sum over A) and refer to Q as the charge 
associated to the corresponding transformation. 

It is well known that the Poisson brackets of the 
conserved charges 7 a define an algebra isomorphic to the 
global continuous symmetry group of the Lagrangian, 
and that the symmetry transformations can be obtained 
in terms of the Poisson brackets of the dynamical vari¬ 
ables {Za) with 7 a. Instead of checking the validity of 
this assertion, later in this section we will use this fact to 
obtain the generators and, from these, the hamiltonian 
for a free particle. 

This is as far as we can get on general grounds. To go fur¬ 
ther we need to specify the symplectic structure. We do 
this taking into account facts known from the mathemat¬ 
ical literature regarding the Galilei group , which we 
assume is the symmetry group of the non-relativistic me¬ 
chanics. It is known that in 3 -I- 1 dimensions, the Galilei 
group accepts only one central charge - to be identified 
with the mass of the particle - and therefore it is not 
possible to introduce NG in 3-1-1 dimensions, at least not 
consistently with the Galilei group. In 2 -f 1 dimensions 
the Galilei algebra accepts three central extensions: the 
mass, the parameter associated to NC and one more that 
we ignore on physical grounds (we are not interested in 
such an extension). Note in particular that consistency 


with Galilei group demands the vanishing of the Pois¬ 
son Bracket among the momenta^Thus, we restraint our 
analysis to 2 -|- 1 dimensions, and assume the symplec¬ 
tic structure given by Eq. 0 . where now 0 = 0 and we 
further assume that: 


dij — dCij , 


(18) 


with 0 a constant parameter, which clearly characterize 
the non-commutativity. The Poisson brackets (PB) and 
the equations of motion read respectively: 


’ ^ dqi dpi dpi dqi dqi dqj ’ 


dH ^ dH 


(19) 


( 20 ) 


Pi =- 


dH 
dqi ’ 


Let us consider now a system invariant under spatial in¬ 
stantaneous transformations: 


St = 0, Sqt = bi, Spi = 0, 


( 21 ) 


where bi (i = 1,2) are arbitrary infinitesimal parameters. 
In order to identify the generators of this transformations 
we write: 

Sqi — {qii Qtrans\ — {Qz; ) (^^) 

Spi — {Pii Qtrans\ — {Pz? — O' 

Thus, using m we conclude that the space translation 
generators satisfy the equations: 


, a. 


dpi 


4 ” dcik ^ — Sij , 


dqk 

9lj 

dqk 


(23) 


= 0 , 


and then, up to a constant, the space translation genera¬ 
tors are the momentum components 7 i = pi. In a similar 
way we can treat the boost transformations: 


St = 0, Sqi = Vit, Spi = mvi- 


(24) 


where Vi {i = 1, 2) are the infinitesimal parameters asso¬ 
ciated to boosts. In this case the transformation of the 
position and momentum are given by: 

— {Qi-)Qboost} — 

Spi = {Pi, Qboost} = {Pz, -kjVj} = mv^. 


^ At this point it is not clear the relation between the Galilei alge¬ 
bra and the symplectic structure. The connection between these 
structures is seen when the generators of the Galilei group are ex¬ 
pressed in terms of the dynamical variables of a physical system, 
see below. 
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Using ini and (El Eq.El reduces to: 


dkj dkj 

"o 

opt oqk 

dkj _ 

dqk " ^ 

therefore the boosts generators h are: 

ki = mqi - Pit + mOeijPj, 


(26) 


(27) 


Likewise, we obtain for the angular momentum, the gen¬ 
erator of spatial rotations: 


J = iijqzPj + -jPkPk- 


(28) 


We can now determine the Hamiltonian of a system pos¬ 
sessing all the symmetries of the Galilean group, which 
we take as the definition of a free particle. To this end 
we consider the time derivative of an arbitrary function 
of the phase space variables: 


dt dt ^ dz°‘ ^ 


(29) 


which by using the equations of motion 0, is cast in the 
form: 


dt 


dj_ 

dt 




Applying this relation to Noether’s charges {Q = 0)we 
obtain the following relations: 


where t is a constant parameter. Therefore, the gener¬ 
ator of time translations is simply 

7 = -H - A. (33) 

Therefore, for every system described by a time- 
independent Hamiltonian, the Hamiltonian itself is a con¬ 
stant of motion. 


For completeness we quote the full Galilei algebra, which 
is obtained through the Poisson brackets and the expres¬ 
sions for the generators previously derived: 


{p^,H} = 0, 
{J,H} = 0, 
{kj,H}=pj, 

[J,ki} = eijkj, 


{Pi,Pj} = 0, 

{JjPi} ~ 

{kj,pi} = mSji, 

{ki,kj} = -m?9eij. 


(34) 


The last two relations show the appearance of the cen¬ 
tral extensions, the mass m and the NC parameter 0. 
Thus, the symplectic structure used to describe 
the mechanics follows from the symmetry of the 
space-time described by the Galilei algebra. 


We now consider the examples at the classical level |^. 
To start let us consider the free particle. The Hamilto¬ 
nian Eq. m leads to the equations: 

Pi . „ 

qi = —, Pi = 0. 

m 

The solution is the the conventional one: 


{H,H} = 0, {p,,H} = 0, (30) 

UH} = 0, {h,H}=p,. 

The Poisson brackets Eq. EU, after using (113, EHi and 
(ESI . amounts to the set of simultaneous equations: 

^ = 0 ^ ^ K 

dqi ’ dpi m' 

Thus, we conclude that the Hamiltonian for a free parti¬ 
cle of mass m, in a non-commutative plane, is given by: 


PiPi 

2m 


(31) 


This is a complicated form of deriving the conventional 
Hamiltonian for a free particle, in this way however we 
are certain that the description is consistent with the 
Galilean relativity. Furthermore, it is important to re¬ 
mark that for a explicitly time-independent Hamiltonian 
H — H{z), the following transformation always define a 
global invariance of the system 


Sz‘^=-^St, 

dt 


qi = diot + qio, Pi = mqio, 

with qi{0) = qio and gi(0) = qio given initial conditions. 
A more involved problem is the system, to which, for 
obvious reasons we will refer as the harmonic oscillator, 
characterized by the Hamiltonian,: 

1 muj'^ 

= -^PiP^ + —^qiqi- ( 35 ) 

In this case the relations arising from the Hamiltonian 
formalism can be combined to produce the following 
equations for x = qi and y = q 2 . 

X = —uj'^x + mOuj'^y (36) 

y = —uj'^y — mOuP'x. 

Denoting the initial conditions by a:(0) = Xq, i;(0) = Xq, 
2/(0) = 2/0 y 2/(0) = 2 / 0 , the solution to the harmonic 
oscillator is given by; 

x(t) = Ti{t)xo + T 2 (t)io+ T^{t) {mOuj '^xq + 2yQ) 

- Ti{t) {mOxo + 2yo) (37) 

y{t) = Ti{t)yo + T 2 {t)yo+ T3{t) [mOuj'^yo - 2xo) 

- Ti{t) {mOyo - 2xo) 


6t = T, 


(32) 



6 


where we have introduced the following definitions: 

cos((/)t) + cos(xi) „ (/'sin(xt) - xsin(^t) 

—- 7 ^ -5 J4(cj — —-— /. , „ 

2 2a;V4 + m2a;^0^ 


^ </'sin(x^) + XsH(t>t) rr. U-, 

= -2^:;^-’ 


cos{xt) — cos{(j)t) 
2a;V4 + rri^uj'^O'^ 


4),X— Y 4" 24= —mOuj^ ^4 + im?9‘^uj‘^, 

As expected, the solution to the NC case, Eq. (EH) re¬ 
duce to the conventional harmonic oscillator in the 0 —> 0 
limit. In order to asses the effect of NC on this system, 
it is convenient to compare the solution to the commuta¬ 
tive case. This is easier to do in the mduj'^ << uj, where 
one shows that: 

x + iy- (38) 

and so we conclude that the effect of NC amounts to 
rotate, in the x — y plane, the commutative solutions 
with angular velocity mQuP'/2. 

To conclude this section we consider the symmetry group 
of the harmonic oscillator. Since the conventional and 
NC cases are different, and can not be obtained one from 
the other (see below), for the sake of clarity we treat 
both of them. The procedure we follow consist in pa¬ 
rameterizing the conserved quantities and then to obtain 
the parameters through the use of the Jacobi identity. In 
fact, if the conserved quantities do not depend explicitly 
of time, then the equation of motion implies: 


Ea. (l42ll to obtain the S'?. The matrix elements 

are constant, unknown parameters, to be determined by 
the relations (13911 y (14011 . After a lengthly calculation one 
concludes that the relations give rise to the following four 
linearly independent matrices : 




(43) 


from which we obtain the conserved quantities: 


S° = 


S? 

Si 


^ {VxVy + , 

1 [2 2 | 22/'2 
^[p^-p^ + mu: {y 

y = ^ i^Pv - yPx) ■ 



(44) 


It is clear that the four quantities are not independent, in 
fact it is easy to show that the following relation holds: 


{iJ,SJ = 0. (39) 

On the other hand, using the Jacobi identity we obtain: 
{H, {S„ •}} - {S„ {H, ■}} = {{H, Si} , •} = 0, (40) 


(^?)' + (^2°)' + (53°)' = ^. (45) 

Moreover one verifies that the S° {i = 1,..., 3) fulfill the 
following algebra: 


Let us begin with the commutative case. We will use 
a subscript or superscript g to avoid confusion with the 
analogous problem in the NC case. The PB of the Hamil¬ 
tonian with an arbitrary function is: 


{4?r}o = - —y— — (41) 

ax m ay opx apy 


We parameterize the PB of the constant of motion 
with an arbitrary function in the following form: 




a 


(42) 


By using this parametrization we assume the conserved 
quantities are bilinear in the phase space variables. Fur¬ 
thermore notice that we have to invert and integrate 


(46) 

Thus we have re-derived the well known result that SU (2) 
is the symmetry group of the two dimensional commuta¬ 
tive harmonic oscillator. 


We now consider the symmetry group for the NC har¬ 
monic oscillator . Instead of Eq. GD) we have to use: 


muX^x 


dpx 


■ muj y 


dpy 


We assume again Eq. (1^ for the PB of the conserved 
quantities. In this case the constraints due to Eqs. (EHl 
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EOli lead only two linearly independent matrices a(j): 


/ 0 —1 0 \ 

1 Om'^Lu^ 0 0 —1 I 

m?u!^ 0 0 0 I ’ 

V 0 0 0 / 

/ 0 1 0 0 \ 

_ I -1 0 0 0 I 
“ 0 0 0 1 ’ 

Vo 0-10/ 


(48) 


so that for the NC oscillator the constants of motion are: 


= H, (49) 

Q 

Si= J = xpy- ypoo + ^{pI+ pI) ■ 

We conclude that SU{2) is not a symmetry of the 
harmonic oscillator on the NC plane The fact 

that diverse conclusions regarding the symmetry of the 
NC armonic oscillator are reached by different authors 
[Sill El is a manifestation of the ambiguities in defin¬ 
ing the theory. A final comment is in order. Comparing 
Eqs. we observe that in the 0^0 limit, the sym¬ 

metry of the commutative harmonic oscillator is not re¬ 
covered. Operationally, it is clear that in the 9^0 limit, 
Eos. 1481149II reproduce only two of the three independent 
conserved quantities. 


classical fundamental degrees of freedom must be Hermi- 
tian and self-adjoint respect to the internal product, and 
second the internal product must be Hermitian, i.e. 

{^\(j))* = {(j)\ilj ), (51) 

the * stands for complex conjugation and |'0), \(j)) are 
arbitrary quantum states. 

Thus, the description of the dynamics of a quantum sys¬ 
tem use the same structures as the classical one, namely 
the Hamiltonian and the symplectic structure. The time 
evolution of the system is given by the Schrodinger equa¬ 
tion: 


Hl-iplt) =ih^\'4’-,t) , (52) 

where H is the quantum Hamiltonian (Hermitian, self- 
adjoint, obtained from the classical time translation gen¬ 
erator) and lip'it) the state of the system at time t. 

From now on we use the notation x = qi,y = q 2 ,Px = Pi 
and Py = P 2 - The commutation relations are then: 

{ [x,x] = [y,y\ = [x,py\ = [y,Px] = 0, 

\Px,Px] = [Py,f>v\ = [Px,Py\ = 0 (53) 

[x,y\ = ihd, [x,px] = [y,Py] = ih. 


III. QUANTUM MECHANICS 

In this section we discuss the problem of quantization as¬ 
suming that a Hamiltonian description - as presented in 
the previous section - is available for the corresponding 
classical system. The quantization proceeds through the 
correspondence principle or canonical quantization 
which associates to the classical phase space a quantum 
space of states described in terms of a Hilbert space. To 
the fundamental degrees of freedom the principle asso¬ 
ciates linear operators z °‘, which act on the Hilbert space. 
The commutation relations of the quantum operators are 
obtained multiplying the classical PB hy ih^. 


{0“, z^} = ^ [i“, = ihrf^, (50) 

Very much as the PB define a geometric structure on 
phase space, the Hilbert space possess an algebraic struc¬ 
ture in the sense that it provides a linea r rep resentation 
of the algebra of the quantum operators [^ . 

The internal product of the Hilbert space must satisfy 
two requirements. First the operators associated to the 


® The h enters naturally in the quantization procedure, just on 
dimensional grounds one requires it. In the NC case the situation 
is different, in fact one should keep in mind the appropriated 
units of 9, which are not the same that at the classical level. 


This is the NC version of the Heisenberg algebra. In this 
section we will develop an abstract representation theory 
of this algebra. To this end we need two basic postulates: 

1. There exists a base \x,Py) that diagonalizes simul¬ 
taneously both the position operator x and the mo¬ 
mentum Py, whose domain of eigenvalues coincides 
with all possible values of the coordinates x and 
Py, which parameterize an arbitrary connected and 
differentiable manifold M. 

X\X,Py) = x\X,Py) , Py\x,Py) =Py\x,Py) , , , 

X, Py G M. 

2. There exists an internal product , positive def¬ 
inite and Hermitian for which the operators x, y, 
Px and Py are self-adjoint. 


With the rules at hand, the application of the postulates 
in the evaluation of the matrix elements {x,py |a;| x',Py') 
and {x,py \py\x',py) imply: 

{x-x'){x,py\x',py) = 0, (55) 

(.Pv-p'y){x,Pv\x',Py) = 0 . 

The general solution to these equations is: 


{x,Py\x',Py) 


S{x - x')S{Py - p'y) 
\/9{.X,Py) 


(56) 









where g{x,py) is a positive definite, arbitrary function 
which is a a priori related to the normalization of the 
eigenbasis. As a consequence of this result the spectral 
decomposition of the identity operator in the \x,py) base 
is: 


1= dxdpy\/g{x,Py)\x,Py) {x,Py\. (57) 

JM 


this completeness relation allow us to construct the wave 
function: 


8 

Evaluation of the matrix elements of the commutators 
{x,y\ and [y,py\, using the postulates and Eg. (15811 . lead 
to the relations: 


{x,py \ [x,y\\x',p'y) =ihe 


5{x-x')5{j}y -py) 


\/g{x,Py) 

= {x - x') {x,Py\y\x' ,p'y) , 


(61) 


V’(a:,py) = (x,py|V’), (58) 

so that for any state \x,py) belonging to the representa¬ 
tion of the algebra lO) 


{x,py\[y,py]\x',p'y) = ih 


d{x - x')d{py - p') 


yjg{x,py) 

= -(Py-p'y) {x,Py \y\ x\Py) . 


(62) 


= / dxdpyJg{x,Py)i:{x,Py)\x,Py), (59) 

Jm ^ 

= dxdpyj g{x,Py)'4}*{x,Py){x,Py\. 

I A/T V 


In particular, the internal product of two arbitrary states 
\tjj) and |(/)), is expressed in terms of the wave functions 
ip{x,py) and (l){x,py): 


{iP\(j))= dxdpyjg{x,py)i^*{x,py)(l){x,py). (60) 
JM ^ 


Thus, matrix element of the y operator can be parame¬ 
terized as follows: 


{x,py \y\x',py) 

= f-0- + —) 

\/9{x,Py) V dx dpy) 

[A{x,Py)+iB{x,Py)\ 


S(x — x')S{py 

x')6{Py-Py). 


(63) 

Py) 


where A{x,py) and B(x,py) are two real, arbitrary functions defined over M. Further constraints on A and B can 
arise from the hermiticity of the y operator, namely: 


—ih 


V9{x,Py) 

ih 


d 


d 




[A{x,Py) - iB{x,Py)\ 

V9ix,Py) 


5{x-x')5{py - p ) 


/^W , [^).x',p'y) + iB{x',p')\ 

/ + T;ir]^ix-x )5{py - Py) + - - 5{x - X )S{py 

^9{x\p'y) V ^gix',p'y) 


-Py)- 


In order to solve this equation we assume the existence of a continuous distribution T{x',Py) on M, we multiply both 
sides of the last equation by T{x',Py) and integrate over the domains of x' and Py 


IM 


\ (-«|- + A) 5(. - - Pi) 


+Was!LlMS$^s{p - r)s{p, - P-) 


\/9{x,Py) 


= / T{x',p' )dx'dp' 


ih 


' M 


- ^')^^Py - Py) 


^J9ix',p'y) V dp'y 


+ 


[A{x' ,p'y) + iB{x' ,p'y)] 

\Jg{x',p'y) 


d{x - x')6{py -p'y) ( , 


(64) 


simplifying this expression and considering that T(x,py) is an arbitrary function, B[x,py) turns out to be: 

^ tV9ix,Py) f ^ ^ 

B[X^Py) - h I 0 + 1 - - 

2 \ dx dpyj y'g{x,Py) 


2 


( 65 ) 
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therefore the matrix elements of y can be expressed as: 


{x,py \y\x',p'y) 


g^/'^{x,py) V dx 


d ^ 5{x - x')5{py - p'y) ^ A{x,py) 
9py) g^/‘^{x,Py) y/g{x,Py) 


x')S(Py 


( 66 ) 


The same approach can be applied to the operator px- 
To this end we consider the matrix elements of the com¬ 
mutators [x,px\ and [px,Py\ to conclude that: 


The explicit evaluation of these relations yields 

p{x,Pv\y\x',p'y)^ = 


-ih d /S{x - x')S{py - p'y)\ 

[x,Py 

g^jiX{x,Py)~^-j-x{x,Pv) 

g^/‘^{x,py) dx V g^/'^ix^Py) J 

ih 

f d\ 

(67) g^'''^(x,py) V ''dx ' dpyj 

-^^(a: X )6{py p ), 

Mx,Py) + [-0-^ + 

Vgi^^Py) 


Vgix^Py) 


x\Py 


Fy) 


g^/^{x,Py) 

-5{x - x')5{py-p'y). 


where C{x,py) is another arbitrary real function defined 
over M. 


An important consequence of the NC version of the 
Heisenberg algebra (151111 arises from the evaluation of 

{x,Py\[y,Px]\x',Py) =0. (68) 

The explicit calculation leads to the compatibility restric¬ 
tion among the functions A(x,py) and C{x,py) 

If xix,Py) is a scalar function defined over M, the trans¬ 
formation 


d 

C{x,Py) - > C{x,Py) - —X{X,Py) (70) 

A{X,Py) -> A{x,Py) + X(x,Py), 

leaves the condition (1^ invariant. This strongly sug¬ 
gests that the functions A(x,py) and C{x,py) are asso¬ 
ciated to the phase definition of the \x,py) basis. Under 
a local U(l) gauge transformation the states transform 
according to: 


and 


p{x,py \px\x',p'y)^ 


= {X,Py 

—ih 


;¥XU,Pt/)p^g--SXU,Py) 


x',Py 


d f 5{x - x')5{py - p' ) 


g'^l^{x,py) dx \ g^A{x,Py) 


+ 


[C{x,Py) - -§^x{x,Py)\ 

V9ix,Py) 


S{x- x')S{py - p'y), 


Thus, under phase transformations of the states, the 
functions A(x,py) and C{x,py) behave according to (17011 . 
Therefore, the configuration space representations of the 
Heisenberg algebra over the manifold M are character¬ 
ized, on one hand, by the function g{x,py), and on the 
other by a flat U(l) bundle defined by the fields A{x,Py) 
and C(x,py). However, since arbitrary local gauge trans¬ 
formations within the U(l) bundle correspond to arbi¬ 
trary local phase redefinitions of the states \x,py), and 
thus relate representations of the Heisenberg algebra 
which are unitarily equivalent, it is clear that all inequiv¬ 
alent representations of the Heisenberg algebra over a 
manifold M are classified in terms of the topologically 
distinct flat U(l) bundles over that manifold, i.e. the 
equivalence classes under local gauge transformations of 
U(l) gauge fields of vanishing field strength over M [^ . 


\x,Py)p = e \x,Py) , (71) 

the matrix elements in the transformed basis are related 
to the original ones according to the following equations: 


.{x,py \y\x',p'y)^ = (x,py 


^^x(x,Pv)y^-T-x(x,Py) 


x',P'y) , 


,(x,Py \Px\x',p'y)^ = (x,Py 


TiX{x,Py)p^^-ix{x,Py) 


X',p'y 


In the case of a simply connected manifold, every holon- 
omy is contractible to the identity. Then, the gauge 
freedom of \x,py) can be used to remove both the 
A{x,py) and C(x,py) fields through the adequate choice 
of x(x,Py). Over a simply connected manifold, the NC 
version of Heisenberg algebra admits only the represen¬ 
tation in which globally A(x,py) = 0 and C(x,py) = 0. 
When the base manifold M is not simply connected 
and therefore possess topological obstructions that pre¬ 
vent some cycles to be contracted, rendering non trivial 
holonomies around them, it is not possible to completely 
remove both A{x,py) and C(x,py). 
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Let us now turn our attention to the wave functions 
'ijj{x,Py) = {x,Py\'ip). Given the parametrization of 
{x,py\y\x\py) and {x,Py\pa:\x',Py), we can use the 
spectral decomposition of the identity operator in order 
to obtain the representation of y and p^ as differential 
operators 


{X,py \y\ ill) = A{x,Py)lll{x,Py) 


(72) 


ih 


g'^/^{x,Py) \ dx 




A. 

dpy 


9^^'^{x,Py)lpix,Py) 


{x,Py IPxIiIj) 


—ih d 
g^/‘^{x,py) dx 


C{x,Py)lil{X,Py) + 
g^/^{X,Py)lll{x,Py) 


(73) 


The last two expressions are the general representation 
of the operators which takes into account: the NC, the 
measure functions g^^^{x,py) related to the inherent ge¬ 
ometry of the base manifold M and to the overall normal¬ 
ization of the \x,Py) basis as well as the possibility of non 
vanishing fields A(x,py) and C(x,py) which could arise 
from the topological properties of M and their possible 
obstruction. 


The properties of the wave function iji{x,py) required so 
that the operators y and px are both Hermitian and self- 
adjoint are 



and 




X,Py) \ {x,Py\lll) I 


2 


= 0 
(74) 




X,Py) \ {X,py\lll) I 


2 


= 0 . 


(75) 


{x,Py\y,Px)- These wave functions are determined by 
the following set of differential equations: 


-ih^ + C{x,Py) 


= Px 


g^/^{x,Py) {x,Py\y,Px) 
g^/‘^[x,py) {x,py\y,px) 


(78) 





+ A{x,Py) 


= y 


g^^'^{x,Py) {x,Py\y,Px) 


g^^'^{x,Py) {x,Py\y,Px) 

(79) 


Since the first order differential equations require only 
one an integration constant, namely the wave function 
(^x^,Py\y,Px), associated to a specific point on the mani¬ 
fold M of coordinates (a;°,p°). Then, any other point of 
coordinates {x,Py) can be reached from {x^,Py) through 
an oriented path P[(x°,p°) i—> {x,py)] running from 
{x^,Py) to (x,py). The solution to (171^1 is of the form 


g^^^{x,py) {x,py\y,px) = g^^^{x°,p°) {x°,p°\y,px) 
X n \p[{x°,p^°^) I —> {x,py)] 


X e 


j- [(x-X°)pa,-(py -Pl)y + 0{py-p^y'> )PX 


(80) 


represents an ordered 


where \j^[{x°,p^'^) i — > ix,Py)] 

holonomy along the path P[{x^°\p^'^) i—s- [x,py)] as 
shown by the following formula: 


^ [P[{x°,pI) '-> [x.Py)]] = 


exp ■ 


+e 


'P{Py-p°) 


' P{x — x'^) 

dpyC{x,Py) 


dxC{x,Py) — 


IP{Py-p°) 


(81) 

dpyA{x,Py) 


Instead of using the basis \x,py), it is possible to work 
with sates \y,Px) which diagonalizes simultaneously the 
position operator y and the momentum component Px ■ 


The normalization condition of the wave function 
{y,Px\y' ,p'x) also requires that 


y\y,Px) = y\y,Px), Px\y,Px) = Px\y,Px), y,Px & D{y,px). 

(76) 

Where D{y,px) stands for the range of spectral values 
of y and px- By analogy with the normalization of the 
\x,py) eigenbasis 15611 . the normalization of \y,Px) is pa¬ 
rameterized according to the relation 


g^P{x°,pl){x°,pl\y,Px) 

so that necessarily 

g^/\x\p^y){x^^\pl\y,Px 


1 

(27rfi)^ yjh{y,pxy 


(82) 


givix°,Py,y,Px) 

{2Trn)h^/^{y,Px)' 


(83) 


{y:Px\y\Px) 


Sjy - y')HPx - Px) 

VHy,Px) 


(77) 


where h{y,px) is again an arbitrary positive definite func¬ 
tion defined over D{y,px). All the results are very similar 
to the ones obtained in the \x,Py) basis, so instead of re¬ 
peating the arguments we discuss the quantities relevant 
to the change among those basis, i.e. the wave functions 


where p{x^,Py,y,Px) is a specific real function indepen¬ 
dent of X and px- Then, the wave functions {x,py\y,px) 
are given by 


{x,Py\y,Px) = 


P^(x°,pl,y,px)^ [P[(x0,p0) ^ (x.Py)]] 

{2'Kh)g^/‘^{x,Py)h^/^{y,Px) 
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The specific choice of ip{x^^\p^,y,px), such that 

^iv{x°,pl,y,px)^-i[x°^px+pi°^y-eplpx] _ ^ (-gg^ 


simplifies the wave function representation of NC Heisen¬ 
berg algebra: 


, , ^ H[P[(x°,p°)i —> {x,py)]]e^^=^P--Pyy+^PyP-'^ 

{2TTh)g^/^{x,py)h^/‘^{y,Px) 

( 86 ) 

This is the NC generalization of the customary wave func¬ 
tion (0 = 0) that arises in conventional QM, and that 
coincides with the Fourier Transform kernel. 

Another admissible basis compatible with the commuta¬ 
tion relations (1^ is Ip) = \px,Py), which diagonalizes 
simultaneously both components of the momentum op¬ 
erator: 

Px \p) = Px \p) , Py Ip) = Py Ip) ■ Px, Py S D{ p). (87) 

Where D{p) is the range of spectral eigenvalues of Px 
and Py. The normalization of|p) can be parameterized 
according to 


This equation is invariant under the following transfor¬ 
mation: 


G.(p) 


Gi(p) -I- 


d^p) 

dpi 


(93) 


with ^(p) as a scalar local function defined on D{p). It is 
important to notice that, in contradistinction to the func¬ 
tions A{x,py) and C{x,py) Eas. ltibIbTII . in the present 
case Gx{p) = 0 and Gy{p) = 0 cannot be simultaneously 
satisfied. For this reason, it is convenient to use a slight 
different parametrization: 


{p\x\p') = 


-k 


{p\y\v') = 




ih f d \ 

d{Px - P'x)d{Pv - p'y) 

l^^^{Px,Py) \dpx) 
(^-lpy + Fx{Px,Py) 

l^^^{Px,Py) 

\ HPx - Px)d{Pv - Py) 

J \/l{Px,Py) 

(94) 

ih f ^ \ 

HPx - Px)S{Py - Py) 

l^^^{Px,Py) \dpy) 

(^^Px+ Fy{px,Py)'^ 

l^l\Px,Py) 

KPx - p'x)d{Py - P’y) 

VliPx^Py) 

(95) 


(p|p') = 


S{Px - p'x)5{py - p'y) 

\/l{Px,Py) 


( 88 ) 


with jiPx-iPy) as a new arbitrary positive definite func¬ 
tion, defined on the D(p) domain. This choice implies 
that the spectral decomposition of the unit operator, in 
terms of the momentum eigenbasis is of the form 


With this convention, Fx(px,Py) and Fy{px,Py) are also 
the components of a vector field defined on D{p) that 
satisfy the condition 


dFyjp) _ dFx{p) ^ ^ 
dpx dpy 

which is invariant under the transformations: 


(96) 


1 = 


lD{p) 


dpxdpy\h{px,Py) |p) (p| . 


(89) 


Fiip) 


FM + 


dm 

dpi 


(97) 


A procedure similar to the one used in the |x,py), leads 
to the following matrix elements: 


{p\x\p) = 


ih 


7^/‘^(Px,Py) \dp: 


d \ 5{px - Px)5{Py - p'y) 


l^l\Px,Py) 


(plyip') = 


, ^ , S{px - Px)SiPy - Py) 

+ Gx(Px,Py) - . . . -(90) 

V'r{Px,Py) 

ih f d \ 5{px - p'x)d{jPy - p'y) 


l^/'^iPx.Py) \dpyj 7^/'‘(Px,Py) 

^S{px - p'x)S{Py - P'y) 

■Gy{Px,Py) - - —■ (91) 

yJlKPx.Py) 


In these equations, Gx{px,Py) and Gy{px,Py) are the 
components of a vector field defined on D(jp), which by 
virtue of the calculation of the matrix elements of the 
relation [x, y\ = ih9, satisfy the compatibility relation 


dGy{p) _ dGx[p) 
dpx dpy 


The p eigenfunction representation of the operators x 
and y are: 

^ 1^1 (^) (98) 

+ ^\py + Fx{Px,Py)'^ '0(p), 

^ 1^1 (^) (99) 

+ (^Px+ Fy{Px,Py)j 1p{p). 


The wave function (x,py|p') can be determined through 
the following set of differential equations 

(Py-Py) {x,Py\Px \p') = 0, 


( 100 ) 
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+ C{x,Py) 


{x,Py\p')^ 

= PW^^(.X,Py){X,Py\p') (101) 




{x,Py\p')) 


= a;7^/'‘(p ) {x,py\p'). 

The normalized solution to these equations is 


( 102 ) 


{X,Py\p) = - 


X 11 


\/2TrHg^^‘^{x,Py)j^/'^{p') 


P[iPx°\pT'‘) '—' ip'x,Py)] 


P[{x°,P^y'’) '-> {X,Py)] 


(103) 


where 12 and ^ are holonomies along the ordered paths 
connecting the fixed points (a;°), p ^'^) and {px^'^ , Py ^^) to 
{x,py) and {p'^^p'y), respectively. 11 was defined before 
in m, and is given by: 


P[ipT\Py°^) '-> ip'x^Py)] 

-i . 40 )) ^P'^+^P(P^ -p(«») 


(104) 


with (l){x^ ,py,px^\px^) as a constant phase. Again, in 
the case of a simply connected base manifold, the gauge 
freedom of the eigenstates \p) can be used to remove com¬ 
pletely the vector field T)(p) through the correct choice 
of the gauge transformation. 


The wave function {x,py\p') can be constructed from 
{x,py\y,px), given by (|8(ill . and from the spectral decom¬ 
position of the unity operator: 


{y,Px\p') = / dxdpyjg{x,Py) {y,Px\x,Py) {x,Py\p') 

J M ^ 

(105) 

\/2TThh^P{x,Py)^^/^{p') 

P[ipT\pT) ' -^ {p'x,Py)] 


where ^{y,px,pT\px°^) 


is an arbitrary real function. 


To conclude this section we quote the change of basis 
among the different wave functions 4'(a:,py) = {x^py\^), 


'^[y.Px) = {y,Px\'^) and «'(p) = (p|4'): 

^{x,py)= f 

dD{y,p^) 

'^ix,Py)= dp^dPyjj{p'^,p'y) {x,Py\p') 4'(p'), 

J D(^Px ■)Py ) 

'^iy,Px)= / dxdpyJg{x,Py) {y,Px\x,Py)^{x,Py), 

Jm ^ 

'^iy,Px)= dp'^dpyj'){p'^,p'y) {y,Px\p) 4'(p'): 

d D{px,Py) 

^{p')= dydpx^/h{y,px) {p'\y,Px) 'i’iy,Px) 

JD(y,P:n) 


(106) 


4'(P') = / dxdpyjg{x,Py){p'\x,Py)'i>{x,Py), 
JM ^ 


with {x,Py\ y,Px ), {x, Py\p') and {y,Px\p') given by equa¬ 
tions ®, ihM and respectively. 


IV. PATH INTEGRAL IN PHASE SPACE 


In this section, we show the relation between the canoni¬ 
cal quantization formalism discussed in the previous part 
and the path integral representation of quantum ampli¬ 
tudes. We follow the conventional approach as well as a 
previous work devoted to the NC case M 


Given a system with initial configuration i, the probabil¬ 
ity associated to the evolution of this system towards a 
final configuration / is 


Pf 




2 


(107) 


If we choose the eigenbasis \x,py) to label initial and final 
states, the transition amplitude can be written as 

K{Xf,Pyf,tf]Xi,Pyi,U) = {Xf,Pyf\e~'^^P~'^''>^ \Xi,Pyi) ■ 

(108) 


The argument is based on the factorization of the tem¬ 
poral evolution operator in the form: 





Inserting two spectral decompositions of the unity op¬ 
erator (in the \x,py) and \y,Px) basis) between each of 
the N factors and making use of the wave functions 
{xj,Pyj\yk,Pxk) Ea. (l86ll . the Kernel can be written as: 
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K (^X f , Pyf j ; ^2 j Pyi t 


(109) 


N- 


= Jim {xf,pyf\ 1 - -eH \xi,pyi 


N 


N-1 


N-1 


= J™ n / dxjdpyj^ 

N^oo IM 


1=1 


fc=0 


dykdpxkX 


—n 


^ dk {Xk-\-l: Pyk+1 \ yki Pxk) {Pk i Pxk \ 1 \xki Pyk) 


N-1 


N-1 


= n / n 

jv^oo A A Jm 


1=1 


fc=0 


dykdpxkX 


:=n 


\/^ 


ri I >■ {xk+i,Pyk+i)] 


{2'Kh)'^g^/‘^{xk+i,Pyk+i)h''-/‘^{yk,Pxk) 
0-1 P[{xl,P°y^,) < - >{Xk,Pyk)] 


1 - 

n 


g^l^{xk,Pyk)h^/‘^{yk,Pxk) 

y^^Ti[{xk + l-Xk)Pxk-{Pyk + l-Pyk)Vk+S(Pyk + l-Pyk)P!ck\ 


where the appropriate boundary conditions for the func¬ 
tional integrals has been indicated using the notation 

x{ti) = Xi,Py{U) = Pyi,x{tf) = Xf,Py{tf) = Pyf. In tWS 

formal expression, the integration measure over phase 
space is the Liouville measure. One can easily identify 
in the last relation the classical action up to an irrelevant 
surface term 



- PyU + OpyPx 


-H] 




UJapZ^ 


= S[z{t)]+A\l'. 


H{z, t) -I- 


dt 


(113) 


It should be clear that the path integral representation 
can be expressed not only using the \x,py) states but in 
any of the basis we analyzed. 


V. WIGNER FUNCTION IN PHASE SPACE 


where, to avoid lengthly expressions, we introduced the 
following notation: 

9] = 9ixj,Pyj), hj = h{yj,pxj), Dk = D{yk,Pxk) 


_ tf-tj , _ {yi,Pxi\H\Xi,Pyi} _ , * 

e- N . ri,- 

i = 0,1, ',iV-I 


Simplifying this expression, we arrive to 

{Xf,Pyf\U{tf,ti)\Xi,Pyi) =0 ^P[{Xj^,Py^) 1 -^ (Xi^Pyi)]^ 


n 


P[ix%P°yf) '-> iXf,Pyf)] 


(2T^Kfg^l\xf,Pyf)g^l\x,,Pyi) ' 


N-1 


N-1 


X lim 
N^oo 


dXjdpy 


M . 


1=1 


'Diyj,px:j) 


dpjdpx 


iX 


X e 


i v^W-1 .1 i^j + l ^j) .. (Pyj + l Pyj) .. I niPyj + i Pyj) ... u 
nl^j =0 ;-yj+t'-;- Pxj—n-: 


As expected, factors 


In this section we discuss the third independent, and 
complete, description of QM, formally distinct to the con¬ 
ventional operator approach in Hilbert space and to the 
Path integral quantization procedure. This quantiza¬ 
tion framework is based on the Wigner quasi-distribution 
function . The main feature of this formalism is the 
fact that interprets the coordinates of phase space z“ not 
as operators but as c-numbers. 

Wigner’s function can be built from the density matrix. 
The density operator is defined as the weighted sum over 
all possible projectors 

p = ^Wn\(l)n) (II4) 

n 

the \4>n) form a complete set of normalized states and 

'^Wn = l. (115) 


n 


PKx^kll,Pyk+l) ' -^ {Xk+l,Pyk+l)] 


X H 


„(0) .„(0) 


P[ixr,N^)^iXk,Pyk)] 


( 111 ) 


in llOhll cancel out among themselves, except for those 
associated to {xf,pyf) and {xi,pyi). 

Finally, the kernel expressed as Functional Integrals over 
phase space (up to irrelevant normalization factors) is: 


{xf-:Pyf \ \Xi,Pyi) 


( 112 ) 


pXf.Pyf t 

/ VxVyVpxVpye^ N -PvV+sPyPx -h] 

' Xi;pyi 


The matrix elements of this operator, with respect to the 
\x,py) basis, for instance, are 

p{x,Py;x',Py) = {x,Py\(j)n) {(l)nW,Py) (116) 

n 

= '^Wn(t)n{.X,Py)(t>l{x' ,p'y). 

n 

Due to the normalization of the states and to Eq. (trrai : 

Tr{p) = J p{x,py]x,py)dxdpy = 1. (117) 
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For a given operator A, the ensemble expected value where a, 6, c, d, f and g are arbitrary constants. Using 

is defined as the Baker-Campbell-Hausdorff formula 


(i) = Tr[pA) (118) 

= ^Wn / (l>n{x,Py){Acl)*^{x,Py)^ dxdpy. 
n '' 

If Wj = 1 and Wi^j = 0 the system is in a pure state; 
otherwise, the system is in a mixed state. The quantity 
Tr{p^) = Yln'^n ^ 1 (with Tr[p^) = 1 only possible for 
pure states), is called Purity. 


The quasi-distribution Wigner function is defined 
through Wigner-Weyl prescription, which assigns a c- 
number function Aw{z) to each operator A in Hilbert 
space. For the two dimensional under consideration, we 
have explicitly: 


^w(^) = I d^ad^T X (119) 

f[{'riPx+r2Py+S-ix+a2y)^ j_ 


xTr 


where we introduced the notation d^ad^r = 
daida 2 dTidT 2 and 4)(cr, r) = Tip^ + T 2 Py + aix + a 2 y- 
The Wigner function W{z) is defined as: 


W{z) = 


1 


(27r?i)4 
xTr 


d^ad^ 






( 120 ) 


g-j;('riPx+T2Py+d-ix+»2y) 


} 


(2 


X {x',p' 


^ / dx'dp' x 

J [Jm 


Q-j[{'^iPx+r2Py+(xix+a2y)^ 


( 122 ) 

the Wigner function is written as: 


xpe 


[2T:hY 

dx'dp'y {x',p'y \ x 

'A)} ■ 


y\ 

p-^<y2y 


This can be still simplified using the operators as trans¬ 
lation generators: 

{e*‘—■»> 

“ T +f l'*l " + ¥ +"f’”»■?)} ■ 


The change of the integration variables C = Ti/2-|-0cr2/2, 
p = —(72/2 is useful to write Wigner function in a more 
compact way 

W(z) = -4^ f dCdple^^^P’^-^y+^'^P-^x (123) 

TT^h^ J I 

X {x - C,Py - p\p\X + C,Py + p)} . 

If the system is in a pure state, with wave function 
'i>{x,Py-,t), then Wigner function takes the form 

= J (124) 

x'i> {x - C,py - p; t) 4-* {x + C,13y + m t)} ■ 


In the following we will restraint to a cartesian and sim¬ 
ply connected NC phase space, so that g{x,py) = 1, 
h(y,py) = 1, 7 (p) = 1 and we can also remove the func¬ 
tions A{x,py), C{x,py) and Fi{p) by means of a local 
gauge transformations. 

Using the commutation relations Eq. and the rep¬ 
resentation of the operators Eq. ™7li and analogous 
equations not explicitly written, it follows that the op¬ 
erators X, y, Px and py, are generators of translations in 
phase space: 

\x,py) = \x +a,py), ( 121 ) 

e~^^y \x,py) = \x + eb,py - b) , 

e-i^^Py \y,px} = \y + c,px ), 

\y,Px) = \y- dd,px - d ), 

\px,Py) = \Px - f,Py) , 

e-isy \px;,py) = \px,py - g ), 


It is also possible to define Wigner function W(z) starting 
from the \y,Px) basis using 11231 IHHll and the spectral 
decomposition of the unity operator; 

W{z) = -4— [ dudvle-^’~^^-Py^+^^Py'>x (125) 

7T^h^ J I 

X {y - v,Px - u\p\y + v,Px -f m)} . 

Similarly, in terms of the \p) basis the Wigner function 
takes the form: 

W{z) = ^ J dudp{e-^i^'^+PP--^^P^P-Py'^^]x (126) 

X {px - U,Py - p\p\Px + U,Py + p)} . 

In analogy to the commutative case {9 = 0), the main 
features of lU)^) are: 

1. Wigner function W{z) is real 

W{z)* = W{z). 
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2. If integrated over x and Py, W^{z) gives the correct 
marginal probability distribution on y and p^ ■ 

\^{y^Px)\^ = j W^{z)dxdpy. 

Similarly, if integrated over y and Px, the Wigner 
function reproduces the probability distribution on 
X and Py: 


This expression can be interpreted as the proof that 
Wigner function cannot be positive definite over 
phase space. If ip and (p are orthogonal, the last 
integral must vanish, and therefore, if W^(z) is not 
equal to zero in a specific region of phase state, then 
it must take negative values in another sector. 

The time dependence on Wigner function follows from: 


Finally, in order to obtain the marginal probability 
distribution on the momentum components, it is 
sufficient to integrate W^},{z) over px and Py-. 

I^(p)l^ = j W^{z)dpxdpy. 

It is important to remark that wave functions dr (p ), 
'^{UtPx) and '^{x,py) are related by means of the 
transformations (lTT7f)ll . 

3. A consequence of the previous feature of Wigner 
function, it is evident that W^i,(z) is normalized 


J W^{z)dydpx 


'd'H (a; - C,Py 


dt 

+ ^ {x-C,Py - 


ih dCdrj r 

- (x + C,Py+ m t) 

d'i/*{x + C,Py + y;ty 

’'■*>- Ft - 


the time dependent Schrodinger equation can be used in 
this expression. In the {x,py) representation the equa¬ 
tion including a potential V(x,y) reads: 


d P 

ih—1p{x,Py-t) = -^%l){x,Py\t) 


fP d‘^'lp{x,Py-,t) 
2m dx^ 


W^{z)dxdydpxdpy = I. 






A 

dpyy_i 


1p{x,Py-,t). 


(127) 


4. Starting from two different density operators pi 
and p 2 , it is possible to construct two different 
Wigner functions Wi{z) and W 2 {z). The oper¬ 
ation Tr{pip 2 ), in terms of Wi{z) and W 2 {z), is 
given by 

Tr[pip2) = (27rft)^ f Wi{z)W2{z)dxdydpxdpy. 


Substituting n?7ii in this relation one shows that, at least 
for potentials which are quadratic in the components of 
the postition operator, the equation reduces to: 

^ Wiz,t) = {H,Wiz)}, 

(128) 


Thus, if Awiz) is a Wigner function associated to 
the operator A HTWi : 


Awiz) = 


1 


d(dT] X 

)} 


TT^h? 

X (x-C,Pv-v\A x + C,py+p 


which is nothing but the familiar Liouville equation ap¬ 
plied to the probability distribution function in phase 
space. Even if the Wigner function does not satisfy all 
the requirements of an authentic probability distribution 
function, it is subject to the same mathematical relations 
as a real one. 


then, the ensemble mean value of A is 


= Tr{pA) = ( 27 r/i)^ 


W {z)Aw{z)dxdydpxdpy. 


VI. QM EXAMPLES ON THE NC PLANE 


5. If a system is in the state and a measurement 
that determines that the new state of the system is 
IP) , then the probability to obtain this result from 
the measurement is |(i/'|^)| . In terms of Wigner 
functions, the transition probability can be written 
as 

|('0|(/))|^ = {2Trhf J Wp{z)W^{z)dxdydpxdpy. 


A. Free Particle 


We begin with the simplest problem, namely the free par¬ 
ticle in two NC dimensions. The corresponding Hamil¬ 
tonian is 

pI+pI 

2m 
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The Schrodinger equation in the momentum representa¬ 
tion determines the spectrum of energies to be the con¬ 
tinuum 




E = 


■pI 


2m 


The general solution is a superposition of stationary 
eigenfunctions of the Hamiltonian: 


il){p;t) = / ase K'^*'iljE{p)dE 


(129) 


However, does not provide any information about 

the wave functions ipip)- order to determine these 
functions, we will take into account Schrodinger equation 
in the \x,py) and \y,Px) basis; 


H'ip{x,py) = - 
Hip{y,Px) = - 


h? d'^'ip{x,py) p\ 


2m dx'^ 




d'^ii){y,Px) ^ 
2m dy'^ 2m 


'<Piy,Px)- 


Using separation of variables to solve those equations we 
obtain: 


ip{x,Pv) 

i’{y,Px) 


'/2TTh 




e^ypv 


(131) 


where <p{px) and fpiPy) are undetermined functions. The 
wave function in momentum space is related to these by 
means of the following transforms 

W)= r ^^S{px - p'x)e-ii<-^P>-^eiyPypipx), 

J_oo 27rft 

= pip'x)e^^P'^P'^S{py-py), (132) 

(133) 

= 4>{p'y)e-~^^P'^P'<‘5{px-p'x)- (134) 

Since ^/’(p) is by assumption separable, then 

'0(p) = dijpx - p'x)S{py - p'y). (135) 

The same result is obtained in the conventional QM, 
which is not surprising, for in both cases Px and Py com¬ 
mute with the Hamiltonian and therefore have common 
eigenfunctions. 


B. Harmonic Oscillator 


Next we consider the 2D Isotropic Harmonic Oscillator, 
described by the Hamiltonian: 

H = ^{pI+pD+ \muj^{dp-+f). (136) 

We will work in the momentum representation of the 
wave function. Since the potential does not depends ex¬ 
plicitly on time, the problem reduces to the eigenvalue 
equation: 


Eil’i.P) = -^{pI + Pl)'^iP) 


ih 


e 

dpx 


-.Pv 


ih 


d 0 

dpy 


^Px 


(137) 


ip{p). 


Rearranging terms, the last equation can be written as 
m^ui'^l 


Etp{p) = ( 1 -f 

h^mu) 


■^{pI + Pl)'^ip) (138) 


2(1 




dpi ^ dpi 


tjj{p) 




On the other hand, the Hamiltonian is rotationally 
invariant, that is, it commutes with the quantum version 
of the angular momentum Eq. P|l : 


0 

J = xpy- ypx + ^{pI+ fy) ■ (139) 

At this point it is convenient to recall that according 
to our analysis of the symmetries of the classical NO 
harmonic oscillator, SU(2) is not a symmetry for this 
system. Thus, for the harmonic oscillator, E[ and J have 
common eigenstates. Noticing that the angular momen¬ 
tum operates over the Harmonic oscillator eigenfunctions 
as follows: 

j^^p) = ,n(^^jy--^Px^i.{p). (140) 

the eigenvalue equation for the Hamiltonian H can be 
written in terms of J as: 

Eip{p) = (^1 + ^ ^ Hoiljip)-^0muj'^J^{p), (141) 

where we introduced the operator 


Hn = 


1,0 o. h^muj'^ 




dpi dpi 
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which is nothing but the Hamiltonian of a commutative The Wigner distribution function is, according to 
harmonic oscillator of frequency w: Eq. 


W = 



(142) 


Besides 


H,J 


= 0, one can shown that 


i7o, J 


= 0 


and 


H,H^ 


= 0 hold. Therefore H, Hq and J have 


common eigenstates . The solution to the partial differ¬ 
ential equation is the product of harmonic oscillator 
eigenfunctions: 


i’nun^ip) 


pI+pI 

g 2mh‘U7 

N 


Hjii ( 


Px 

mhzu 


)HuA 


Py ) 
mhzu'^ 


(143) 


where N = Thus, the eigenvalue 

equation for Hq is: 


Holism,n^ip) = hw {m + n2 + 1) 4^ ni,n2 (P), (144) 


Wq^q{z) = 


'J TT^mh^TU 

/ OO 2 2 

dudije-"^ e-^ iPpy-Pv'^)] , 

-OO 

1 {-4i#-^[(-+fp«)'+(y-fp»)']} 


■ (148) 


This function is positive definite in its entire domain. 
Without loss of generality let us consider for simplicity 
the case m = l,h = l,u} = 1 (which can be obtained 
through a scale transformation). The time evolution of 
Wigner function is determined by Eq. The solution 

to this equation can be written operationally as: 


Wo, 0 ( 2 :, t) = exp [it {H, •}] Wo, 0 ( 2 ,0), (149) 

which is nothing but a linear, time dependent transfor¬ 
mation. Then, the complete solution is 


with ni,n 2 positive integers. These functions form a 
complete orthonormal basis. The lineal combination of 
these functions that simultaneously diagonalizes H and 
J is: 


P^'^Py ^+1 i 


i^n,j{p) = 


g 2Tnh'co 

w 


EE 

r—O q—0 


^+3 


^-3 


(145) 


m-i)»«««//„-<,«,(^)h.«(=;) 


\mhw J 


with n G , 2j G Z, subjects to the restriction 
-f < j < t and iV = 2^V7rmtiu7y/(f +j)! (f-j)'- 
Putting all together the eigenvalue spectrum associated 
to eigenfunctions (I145II is summarized in the equations: 


H-ipnjip) = Enj-ipudip) 


-1^2/ .2^2 

Hujxll -\ ---(n -I- 1) — Omuj'^hj 


(146) 

i’njip), 


Wo,o( 2 , t) = exp [it {H, •}] X 


(150) 




exp < -- 


V4-f 


{pI+pD \ X 


X exp ■ 


{x+iPvf+{y-y<py 




C. 2D Solid (Einstein’s model) 


The thermodynamic properties of this solid can be stud¬ 
ied considering a canonical ensemble of N distinguish¬ 
able, non interacting and NC harmonic oscillators [^ . 
For an oscillator of frequency w, the probability of being 
in the n,j state, denoted Wnj-, is 


Wn,j 


Z(T,E,1)’ 


(151) 


Ji’njip) = 2hjll;ri,j{p)- 


To conclude the discussion on the harmonic oscillator, 
we derive the Wigner quasi-distribution function for the 
ground state of the harmonic oscillator, which is de¬ 
scribed by the wave functions 


g 2mhz!7 

V'O.o(p) = y ^ 
yirmfiw 


(147) 


where Z{T, V, 1) is the partition function 

°° E 

Z{T,V,1) =Tr{e-^) = Y^ ^ e"^. (152) 

n=0 j= — n/2 


The calculation of this function proceeds through geo¬ 
metric sums, and renders the following result: 


Z{T, V, 1) 


= 2 


cos 


f hujQ\ f twj'^md 

-TTT^ - cos _ 


\2kT J 


2kT 



( 153 ) 
















































18 


with 


is given by 


0 = y/A + m^Lu^O^. (154) 

At this point, it is necessary to know the distribution 
function of the natural frequencies in the solid. To sim¬ 
plify the calculation we choose Einstein’s approximation, 
and set all frequencies equal LOi = uj. This approach, if 
not the better, provides a clear qualitative behavior of 
the system. 

The partition function of the N oscillators is, then, 

Z{T,V,N) = Z{T,V,l)^. (155) 


dU 


(161) 


V,N 


h^uj'^N 

2fcT2 

(2 + 

—mojOQ 


cos 


fhujQ\ 
\^) ~ 


cos 


2 2/)2 
m UJ “ 


9^) 


cos 


I 2fcT ) 


sm 


/hujQ\ 

) 


sm 


cos 

huj^mO 

2kT 


huj^mB 
2kT 
hw^mO 
2kT 


1 -2 


- 1 


Again, in the high temperature regime one recovers the 
conventional result: 


From this partition function, the derivation of the Free 
Energy is immediate 


Cv = 2kN - 


fi^uj'^N (2 -I- m?uj‘^6^) 
i2fcr2 


(162) 


A{T,V,N) =-kT\n [Z{T,V,N)], 

/ hujQ\ 


(156) 


V2fcr ) 


= NkT\n cos 
+ NkT \n2. 

Thus, the entropy of the system is 


— cos 


huj^mO 

2kT 


Finally in figures ^ and |21 we show the behavior of the 
entropy 5" as a function of T and 0, for fixed values of 
m, UJ and N. In the present model, the entropy varies 
significatively with respect to 9, in particular if 0 ~ 
{muj)~^. 


dT 

hujN 


= -Nk\n 


V,N 


COS 


/ftw0\ 

V2 ^) ~ 


COS 


fvjj^m9 

2kT 


V4 + m2w202 sin (|||) - mw0 sin 


2r[cos(tf)-cos(i^)] 

-7Vfcln2, (157) 

and its internal energy takes the form 

U = A + TS (158) 


hjjjN 


0 sin (Iff) - muj9 sin 


2 T [cos (Iff)-cos (^)] 


VII. SUMMARY AND CONCLUSIONS 


We presented a systematic study of non-commutative 
mechanics starting from the classical formalism and pro¬ 
ceeding through the quantization. We emphasized the 
role played by symmetries, in particular we ensure that 
the NC free particle is consistent with Galilean relativity 
which follows from the well known relation among the 
boost generators and the position operator. The gen¬ 
eral description of NCCM in terms of second class con¬ 
strained system was elaborated for Hamiltonians of the 
type H = T + V where T and V stand for kinetic and 
potential energy of the two dimensional system. 


In the high temperature regime kT ^ hu, the behavior Besides providing a global view of the problem, our 

of U can be deduced expanding in power series of T : manuscript contains new results, in particular: 


UkT:^hui — 2NkT -\- 


h^uj'^N (2 -I- m^uj‘^9‘^) 

mr 


(159) 


We conclude that the conventional energy equipartition 
is also obtained in the NC case. On the other hand, in 
the opposite limit T —> 0, the internal energy reduces to 
the contribution of the minimum energy of each oscillator 


lim U = Nhjj 

T^O 



m?uj'^9'^ 

4 


(160) 


• A formulation that avoids the use of non-canonical 
transformations and/or expansion in the NC pa¬ 
rameter 9. 

• In classical mechanics, analytical solution for the 
free particle and harmonic oscillator problems. We 
also show that SU (2) is not the symmetry group of 
the isotropic harmonic oscillator, and identify the 
generators of the existing symmetry. 

• Quantization is presented in three different frame¬ 
works: Canonical, Path Integral and Wigner Func¬ 
tion. 


Finally, the calorific capacity of this set of NC oscillators 
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• The representations of the Heisenberg Algebra in 
three different basis {{Px,Py)^ i^jPy) {Px^u)), 
including the wave functions permitting the change 
of basis. 

• Non-equivalent representations of the Heisenberg 
algebra characterized by gauge fields that follow 
naturally from the structure of the algebra. Those 
Fields are relevant in the description of multiply 
connected manifolds. 

• Extension of the analysis of the fundamental prop¬ 
erties of the Wigner function in four dimensional 
NC phase space. 

• Analytical solution in QM for the free particle and 
harmonic oscillator without performing the cus¬ 
tomary non-canonical transformation, without us¬ 
ing the structure of SU (2) generators and without 
assuming a priori the existence of a vacuum state. 

• The thermodynamic properties of a 2D NC crystal 
using the 2D Einstein’s solid model. The behavior 
of the entropy as a function of the NC parameter 9 
is reported. 

The proper definition of the system to be treated in 
NC mechanics and the certitude that it is not traded 
along the way are fundamental. Our approach focuses on 
both points, first defining a general consistent framework 
and second avoiding completely the use of non canonical 
transformations since those lead, in general, to a system 
whose properties are completely different from those of 
the original one. 
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FIGURE GAPTION. 

Figure 1. Normalized entropy as a function of the NG 
parameter 9 and temperature as obtained from Eq. ilFfll . 

Figure 2. Normalized entropy as a function of temper¬ 
ature, for fixed values of the NG parameter 9 as obtained 
from Eq. 
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